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SOME THEORETICAL LOW-SPEED SPAN LOADING CHARACTERISTICS OF 
SWEPT WINGS IN ROLL AND SIDESLIP 

By John D. Bibd 


SUMMARY 


SYMBOLS 


The Weissinger method for determining additional span 
loading for incompressible flow is used to find the damping in 
roll, the lateral center of pressure of the rolling load, and the 
span loading coefficients caused by rolling for wing plan forms 
of various aspect ratios, taper ratios, and sweep angles. In 
addition, the applicability of the method to the determination 
of certain other aerodynamic derivatives is investigated, and 
corrections for the first-order effects of compressibility are 
indicated. 

The agreement obtained between experimentally and theo- 
retically determined values for the aerodynamic coefficients 
indicates that the method of Weissinger is well suited to the 
calculation of the additional span loading caused by rolling 
and for the calculation of such resulting aerodynamic deriva- 
tives of wings as do not involve considerations of tip suction. 

INTRODUCTION 

The National Advisory Committee for Aeronautics has 
conducted for some time a program consisting of both 
theoretical and experimental investigations for determining 
the effects of plan form on the various aerodynamic deriva- 
tives of w ings . The efforts in this direction have recently 
been intensified as a result of recourse to large amounts of 
wing sweep for delaying the onset of the drag rise associated 
with Mach number effects. Several methods for computing 
some of the low-speed characteristics of wings having large 
amounts of sweep have been evaluated in reference 1, and 
one of the most promising methods was used in reference 2 
for an analysis of the effects of plan form on the additional 
loading caused by angle of attack. 

In this report, the low-speed effects of plan form on the 
additional loading caused by rolling and on the damping in 
roll are deter min ed. The calculations were made by' the 
method of Weissinger (reference 3) for a variety of plan 
forms having sweep angles, aspect ratios, and taper ratios 
which encompass the probable ranges of these variables. 
Experimental data are given for verification of some of the 
results, and the applicability' of the method to the deter- 
mination of certain other derivatives is investigated. The 
simplification of the Weissinger method resulting from the 
fact that the span loadings under investigation are anti- 
symmetrical about the midspan of the wing is given in the 
appendix along with a calculation form and necessary 
instructions for its use in calculating antisymmetrical lift 
distributions. Constants are given for calculations which 
determine the circulation at 7 and 15 points along the wing 
span. 


Forces and moments are referred to the stability' axes, the 
origin of which is assumed to be at the projection of the 
quarter-chord point of the mean aerodynamic chord of the 
wing on the plane of s ymm etry. The stability-axes sy'stem 
is shown in figure 1. The coefficients and symbols used 
herein are defined as follows : 
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spanwise eenter-of-pressure location of one wing 
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rolling moment 
yawing moment 

angle of attack of wing with respect to axis system 
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mean wing chord 



angle of sideslip, radians 

dihedral angle (measured in plane which includes 
F-axis and is normal to wing surface), radians 
rate of roll, radians per unit time 
dynamic pressure 
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A angle of sweep of quarter-chord line of wing 

„ x x . / Tip chord \ 

X taper ratio 

y distance along F-axis from origin 

y ep distance along F-axis to lateral center of pressure of 
loading under consideration 

do section lift-curve slope for section normal to quarter- 
chord line, per radian 

M 0 Mach number, ratio of velocity of free stream to 
velocity of sound 

RESULTS AND DISCUSSION 

GENERAL REMARKS 

The calculations of span loadings carried out in this report 
were made by the method of Weissinger as described in 
references 1 and 3 and as utilized in reference 2 for an investi- 
gation of the effects of plan form on the additional loading 
causedjby jangle of attack. The span loadings caused by 
rolling, the resulting damping-in-roll parameter C ip , and the 
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Figure 1.— Stability system of axes. Positive values of forces, moments, and angles are 
indicated by arrows. 


lateral center of pressure ij ej , are given for the range of plan 
forms shown in figure 2 along with a check of the validity of 
the method of Weissinger for rolling loads in figures 3 to 6. 
The values of O n JC L and the dihedral effect contributed 
by geometric dihedral (7^/r, as determined from the calcu- 
lated loadings and related information, are presented in 
figures 7 to 9 for a few untapered plan forms. 

The span loadings presented herein and the values of the 
aerodynamic derivatives may be corrected to correspond to 
section lift-curve slopes other than 2iir by multiplication of 
the loading coefficients or derivatives by the ratio of the 
chosen lift-curve slope to 2x (which is inherent in the Weis- 
singer method) as is discussed in reference 1. A somewhat 
more exact way to account for the effect of changes in the 
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Figure 2 . — Range of plan forms for which rolling loads were calculated. 
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slope of the section lift curve on £7i would he to obtain the 
following ratio (reference 4) : 



A + 4 cos A 


(^) ST A+4 cos A 


Multiplication of the Weissinger values of £7^ by this ratio 
would vield the value of £7, for any chosen section-lift-curve 
slope. 

The subsonic effects of compressibility as given by an 
application of the Prandtl-Glauert rule may be determined 
for the characteristics presented herein in the same manner 
as given in reference 2 for the additional load caused by 
angle of attack. The characteristics of a wing at a given 
Mach number are obtained for an equivalent wing in incom- 
pressible flow as defined by 
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■Anuhalent — V 1 — AIq 2 
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The loading coefficients and derivatives 
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and Cip/r are then multiplied by l/Vl— -W- The lateral 
center of pressure 17 ^ needs no multiplication factor. 


COMPARISON OF RESULTS FOR ROLL WITH EXPERIMENT 


Figure 3 gives a comparison of values of C tp as calculated 
by the 7-point method of Weissinger and as determined by 
experiment in the rolling-flow section of the Langley stability 
t unn el. The values indicate very good correlation for most 
of the test results when the section lift-curve slope of the 
experimental data is taken into account. Figure 3 indicates 
that the section lift-curve slope of the experimental data is 
90 to 95 percent of the theoretical value of 2x. The method 
of Weissinger is expected to apply well to the case of rolling 
loads because the control points for the downwash sum- 
mation are grouped near the tip of the wing where the load 
is greatest and is changing most rapidly. Calculation of 
rolling span loadings by the 15-point Weissinger method 
results in no appreciable improvement over the 7-point 
method for the range of plan forms investigated. The 
accuracy of the Weissinger method was found, however, to 
decrease with increasing aspect ratio and sweep. The 
7-point method is used herein to determine the effect of plan 
form on the rolling span loading, the damping in roll, and 
the lateral center of pressure. 


EFEECTS OF PLAN FORM ON LOAD DISTRIBUTION AND DAMPING IN ROLL 

The additional loadings caused by rolling, the damping in 
roll, and the lateral center of pressure of the additional 


loading for wings of various aspect ratios, sweeps, and taper 
ratios as determined by the 7-point method of Weissinger 
are presented in figures 4 to 6. 

For aspect ratios of the order of 6 the lateral centers of 
pressure move outward, and the values of £7j # decrease with 
angle of sweep in a normally accepted manner for the range 
of plan forms investigated. The variation of the loading 
along the span becomes markedly more linear for a greater 
proportion of the span with increased sweep for taper ratios 
of 0.5 and 1.0, this linearity indicating an approach to the 
case of the two-dimensional yawed wing. A few calculations 
for negative angles of sweep at an aspect ratio of 3.5 and a 
taper ratio of l.'O show that the value of for a positive 
angle of sweep is very nearly the same as that for the cor- 
responding negative angle of sweep for this particular plan 
form. The effects of increased taper (decreased taper ratio) 
are to shift the load toward the center of the wing and also 
to decrease the magnitude of the load for a given rolling 
velocity. These effects, of course, result in decreased values 
of C t . For the range of taper ratios considered herein, most 

V 

of the reduction in C t with taper ratio occurs between taper 
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ratios of 0.5 and 0. 
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Figube 3, — Comparison of values of C\^ determined from JTACA tests and by 7-point 
method of "Weissinger for a range of wing plan forms. Lines of perfect correlation are 
shown for two ratios of experimental to theoretical section lift-cnrve slope. 
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Figure 4.— Additional span loading in roll for wings of various tapei' ratios and sweeps as calculated by the 7-point method of Weissinger. 
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Fiqttbe 5.— Values of Cj 
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(a) X=0. 

(b) X=0.5. 

(c) X=1.0. 

for wings of various taper ratios, sweeps, and aspect ratios as escalated by tbe 7-point method of Weissinger. 
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Taper-ratio and sweep variations produce little effect on 
the magnitude or distribution of loading for an aspect ratio 
of 1.5 and, consequently, produce small changes in the values 
of Ci v and the lateral center of pressure t\ cv . Reference d 
indicates a reversal in the effect of sweep on the lift-curve 
slope at low aspect ratios, and sweeps. A similar effect was 
noted herein for C^, but the magnitude of this reversal was 
exceedingly small. 

The result of reference 5, for extremely Iow-aspect-ratio 
triangles, is included in the plots of C* against aspect ratio A. 
(See fig. 5.) This result in conjunction with the material of 
this investigation indicates a logical approach to zero plan- 
form effect at zero aspect ratio. Reference 2 indicates that 
this result is true for the lift-curve slope also. 

COMPARISON OF RESULTS FOR C. WITH EXPERIMENT 

V 

The results of a few calculations of the derivative C„ for a 

V 

taper ratio of 1 .0 are presented in figure 7. These calculations 
involved the use of span loadings determined by the method 
of Weissinger and the concepts of the Prandtl lifting-line 
theory. The additional loading caused by angle of attack 
was determined from reference 2 and the loading caused by 
rolling, from the present report. The procedure, in effect, 
involved finding the direction cosines of the lift vector, at 
various stations along the wing in roll at an angle of attack, on 
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Figure 6.— Lateral center of pressure of additional loading caused by rolling for wings of . 
various taper ratios, sweeps, and aspect ratios as calculated by the 7-point method of 
Weissinger. , 


the assumption that this vector was at all times perpendicular 
to the relative wind and the quarter-chord line of the wing. 
The final force and moment coefficients were determined by 
integrating the components of the vectors across the span. 
This method of calculation which accounts for the effect of 
leading-edge suction appears to be a better approximation 
for the rectangular case than for highly tapered wings. 

The results of these calculations are shown in figure 7 
along with experimental data from reference 6 and an extra- 
polation of the results of reference 7 as presented in refer- 
ence 4. The discrepancy between the results of these calcula- 
tions and the experimental points appears to be attributable 
to the neglect of tip suction. The contribution of tip suction 
to 0 Up is known to increase as the aspect ratio is diminished. 
For the subsonic case, as the aspect ratio of a rectangular wing 
is made smaller, more of the lift is carried at the leading edge. 
For vanishingly small aspect ratio, Ribner's theory (refer- 
ence 5) may be extended to give the value of C n JC L for the 
rectangle. This result is shown in figure 7 by the dash line. 
Good agreement is obtained with the experimental results 
for the rectangle. This theoretical result is purely the effect 
of tip suction, and its. agreement with the experimental 
values indicates that C n ^jC L is largely dependent on this 
phenomenon for Iowaspect ratios. The method of Weissinger, 
which contains no consideration of chord loading, thus 
does not appear to be ideally suited to the calculation of 
the increments of derivatives contributed by tip suction. 

COMPARISON OF RESULTS FOR C tf lr WITH EXPERIMENT 

The results of the application of the 15-point method of 
Weissinger to the calculation of the effect of 100-percent-span 
geometric dihedral on the additional span loading caused by 
sideslip and the derivative C if are given in figures 8 and 9. 
In calculating the loading caused by sideslip, the effect of the 
skewness of the trailing vortices was assumed to be small in 
comparison with other considerations. This assumption 



Fioure 7.— Comparison of experimental and calculated values of Cn 1> IC L . \-l.o. 







Fig cee 9.— Results of calculation of C^f V for wings of various amounts of sweep by the 
15-point method of Weissinger and comparison, with experimental data. X= i.O. 


reduces the problem to one of determining the span loading 
of a wing which has equal and opposite angles of attack on 
the two wing panels. These angles correspond to the 
product of the angle of sideslip and the angle of dihedral for 
the right wing panel and the negative of this product for the 
left wing panel. 

Comparison in figure 9 of the results of the calculation of 
Ci (T by this method with experimental results from refer- 
ences 8 to 13 indicates good agreement in the case of zero 
sweep. Experimental results are given in some cases for 


several angles of attack. It should be noted that this 
calculation predicts appreciably smaller values of C t jY 
than the results of reference 7 and the extrapolation thereof 
(reference 8) shown in figure 9. These results approach 
those of reference 5 for low-aspect-ratio triangles as the 
aspect ratio is diminished. The effect of increased sweep is 
to decrease the value of C t jT, the effect being greater for a 
given incremental sweep at large sweep angles than at small. 
The approximate method of reference 8 which was designed 
to predict changes in C t jY with changes in sweep predicts 
the effect of sweep to be greater than that shown in the 
calculations of the present report. 

The additional loadings caused by a unit dihedral at unit 
sideslip indicate a concentration of the loading much nearer 
the center section than in the rolling case which is in agree- 
ment with the anticipated result. (See fig. 8.) The effect 
of sweep is again to shift the load nearer the tip for the range 
of plan forms investigated. 

The 7-point method of Weissinger predicted slightly 
smaller values of C t f r than the 15-point method. This 
discrepancy is in all likelihood a result of the lack of definition 
of the unit angle of attack in the vicinity of the midspan. 
The Weissinger method groups most of the control points for 
the summation of the downwash near the tip region. For 
calculations of the type discussed in this section and for 
aileron or twist investigations, a uniform spanwise grouping 
of the control points appears to be preferable. 
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CONCLUSIONS 

A theoretical investigation of the low-speed additional 
span loading of wings of various plan forms in roll and certain 
related derivatives indicated the following conclusions: 

1. The method of Weissinger appears to be well suited to 
the calculation of the span loading caused by roll and for the 
calculation of such resulting aerodynamic derivatives of 
wings as do not involve considerations of tip suction. 

2. For aspect ratios of the order of 6 the lateral centers of 
pressure move outward, and the values of C t , the rate of 
change of rolling-moment coefficient with wing-tip helix 
angle, decrease with angle of wing sweep in a normally 
accepted manner. The effects of increased taper are, of 
course, to shift the lateral center of pressure toward the 
center of the wing and to decrease the values of C t . 

P 


3. Taper ratio and sweep variations produce little effect 
on the magnitude or distribution of span loading for wings of 
aspect ratios of 1.5 or lower and, consequently, produce small 
changes in the values of the damping in roll and the lateral 
center of pressure of the rolling load. 

4. The effect of increased sweep is to decrease the value of 
the rate of change of rolling-moment coefficient with angle of 
sideslip per unit dihedral angle, C t jT, the rate of decrease 
being largest at large sweep angles. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Air Force Base, Va., December 22, 19//.8. 


APPENDIX 


PROCEDURE FOR CALCULATING ANTISYMMETRICAL LOADS BY THE METHOD OF WEISSINGER 


GENERAL 

Reference 1 presents a development of the method of 
Weissinger and a procedure for calculating symmetrical 
loads. The calculation of antisymmetrical loads permits, 
in much the some manner as in the case for symmetrical 
loads, an appreciable reduction in labor over that of the 
general solution. These simplifications and the tabular form 
used for the commutations made for this report are presented 
herein for completeness. The use of the tables for calculat- 
ing span loadings is described at length in a section following 
the derivation of the equations. 

DERIVATION OF EQUATIONS 

The method of W eissinger considers the wing replaced by 
a bound vortex located on the quarter-chord line and trailing 
vortices extending downstream to infinity. The distribu- 
tion of circulation is determined so that the vertical com- 
ponent of the induced velocity due to this vortex system at 

the -qC line is equal and opposite to the corresponding com- 
ponent of the incident flow. The required integrations of 
the downwash equations are performed by summations in 
the manner of Multhopp (reference 14). Weissinger derives 
a set of simultaneous equations in terms of the loading 
coefficient O n (reference 3) 


where m is the number of points at which the span-loading 
coefficients are to be determined. The symbol a, refers to 
the geometric angle of attack at station v, and the loading 
coefficient G n is the circulation at station n divided by the 
wing span and free-stream velocity. The coefficients &*„,« 
in this set of equations are given by 

b*r,n—%b V ' n — ar v g, -n (2) 

Also, 

f me \ 

h sm Wi; n-(-i)— q 

fe)J L 2(m+1) J 

except for n=v where 


m+1 

4sin (svi) 


The parameter ar v is the ratio of the wing span to the chord 
at station v. The influence function g v , n is defined by 

_ i r m 

9r ’ n ~2(M+l) LS L ^ v ’ 


— S b*v,nGn = CCy 

71=1 


( 1 ) 


La (r , 0)/ B ,o + L a (r , Af + 1)/« 
2 


o 


(3) 


wherein the integers v aiid n are related to the nondimensional 
spanwise stations i? and ij by 


and 



VTT 

m+1 


;=J L = , 
6/2 


nv 

m+1 


and is the result of an integration over the span of the plan- 
form function La multiplied by the values of/: 

r 2 J2, me nie ... 

^■'"m+1 2 ^ sm ^ 5+1 “ I+I 

where M is the number of stations used in the spanwise 
integration of equation (3). 
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For an antisymmetrical load 

(on-J-l— * 


G m 4-1=0 
2 


and 


CCy — P 

a m+ 1=0 
2 

Tims, the solution of the downwash equation now becomes, 
because of the reduction in the number of unknowns, 


Where is odd, the functions containing and /m+i-K of 
equation (7) cancel, so odd values of /i x may be eliminated 
from further consideration. Further examination of equa- 
tion (4) shows that 

Applying equations (8) and (9) to equation (7) gives 

TO— 1 / 

^ j a) + 

7t=l t 


— ’y',(b* r . n — b*,.m4-l-n)Gn=aL, 

71 = 1 

M - 1 


( 6 ) 


If the summation is made from 1 to equation (3) 

becomes 


J. 

9r ‘ n=z ~2(M+l) 
-LaO, Af+1 



[Aa(*',Ai)t 

-Aa(i-,M+1-ai)]= 






;nS[AA(^)/n,,+ 

«v^/l 

i+tan A(,-«] 

+ (’?— n) 2 ~l 

4- 


L\(v, OI/k.q+T'aCp, 
2 


gives the solution 

yr 

— I %(br,n ») i 

n=i L 


n.Af+l ^ 


(»^±. l )/.3£>+ 


ar r (ri — ij) 


1 - 

[Vr^: +tanA(7J ~’ ) _ 

+ (’7+’7) 2 ^ 

^ (6) 

ar,(v+v) 

l . tan A J 

id (6) into equation (5) 

2 tan A-y/ 

1 

or 

(aF+ v tanA ) ^ 

; — {— 2?7 tan A 

V 



/' iV/4-l\ r , iA(j',0)/ n>O +7 | A(!',M'T-l)/ a ,.vf+I 

Aa -y- j / n ,M±l+ 2 

| S [Aa ( y , n)jm+l~n , M + Aa(» , M+ 1 — Jf+l-jt] 1 + 

r / Af+1\ ^ r Aa(p, 0)./ m +i- n ,o , 

Aa — g— J Hi-i-i 2 ~ T 


Ai(v , M-f l)/rn+l- 


An examin ation of equation (4) shows that for even values 
of Jit 

/m+l-n,M = '/rc.lt (®) 

and for odd values of y x 

3m+l—n,n = Jn,ii 


M + 1 


2{M+l) — ^ a (jo m) + A A (y 


, M-i-1— ft)] j<?n=av 


(10) 


JUT 1 1 

where ^ has one-half its normal value for /i=0 and — ^ — 

Since ^+i-a=— ^(i“) and only positive values of rj are 
used in this calculation form, from reference 1, page 51, 


+ 


For convenience, the parameters 

2 (bv, n b Vt m +.i—n) 

[AaCv, m) + AaC*', M+ 1 — it)] 

and 

M+l 

S — 2/ re ,^ [A A (r, m) T-A A (r,il/+l— m)] 

^=o 

are denoted by 


and 


fn,P 

If 

***, n 


respectively. Equation (10) thus becomes 


TO— I 


S [ C, ’ n+ 2(M-f 1) 9, ’"1 <?n-a ' (n) 



388 


REPORT 969 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


USE OF TABLES 


REFERENCES 


Table I presents the tabular form used for making the 
calculations where m=M=7. Table II presents values for 
certain columns of table I for use when making calculations 
where m=M= 15. The 15-point calculation takes approxi- 
mately four times as long to. complete as the 7-point 
calculation. 

A spanloading is calculated by first selecting the wing sweep, 
aspect ratio, and taper ratio and carrying out the operations 
indicated in table I through column ©. Column © is then 
calculated by the procedure given in the sample calcula- 
tions in table III. The remainder of the columns of table I 
are then calculated as indicated. The values of column ® 
are the coefficients of the loading ordinates G n for v=l, 2, 3 
. and /i=l, 2, 3 as required by equation (11). These values 
determine three simultaneous equations relating the loading 
at three spanwise stations (corresponding to the values of 
n) to the values of the angle of attack at three spanwise 
stations (corresponding to the values of v). For determining 
the loading caused by roll (unit pb/2V), these angle-of-attack 
values are the same as the values of ij given in column ®. 
Table III shows the setup of the simultaneous equations 
from the values of column @. These equations can be most 
easily solved by the method given in reference 15. The 
results are given in the form of loading ordinates G„ which 
refer to the spanwise positions given by 

— ■ y nr 

v ^tj2 =cos m+l 


The loadings at the positions n = 0 and 


m+1 

2 


, corresponding 


to the tip and the root, are zero. 

Calculations made for m=M= 15 are carried out in much 
the same way as described herein for m—M= 7. The 
columns of table II should be substituted for the correspond- 
ing columns of table I before proceeding with the calculation, 
however. This calculation yields seven simultaneous equa- 
tions for solution and a corresponding number of loading 
ordinates. 

The rolling-moment coefficient resulting from known load- 
ing ordinates may be determined for antisymmetrical loads 
by 


m— I 


irA 


2(m+l) 


S G n sin 


2mr 

m+1 


1. Van Dorn, Nicholas H., and DeYoung, John: A Comparison of 

Three Theoretical Methods of Calculating Span Load Distribu- 
tion on Swept Wings. NACA TN 1476, 1947. 

2. DeYoung, John: Theoretical Additional Span Loading Charac- 

teristics of Wings with Arbitrary Sweep, Aspect Ratio, and 
Taper Ratio. NACA TN 1491, 1947. 

3. Weissinger, J.: The Lift Distribution of Swept-Back Wings 

NACA TM 1120, 1947. 

4. Toll, Thomas A., and Queijo, M. J.: Approximate Relations aud 

Charts for Low-Speed Stability Derivatives of Swept Wings. 
NACA TN 1581, 1948. 

5. Ribner, Herbert S.: The Stability Derivatives of Low- Aspect- 

Ratio Triangular Wings at Subsonic and Supersonic Speeds. 
NACA TN 1423, 1947. 

6. Goodman, Alex, and Fisher, Lewis R.: Investigation at Low 

Speeds of the Effect of Aspect Ratio and Sweep on Rolling 
Stability Derivatives of Untapered Wings. NACA Rep. 968, 
1950. 

7. Pearson, Henry A., and Jones, Robert T.: Theoretical Stability 

and Control Characteristics of Wings with Various Amounts of 
Taper and Twist. NACA Rep. 635, 1938. 

8. Queijo, M. J., and Jaquet, Byron M.: Investigation of Effects of 

Geometric Dihedral on Low-Speed Static Stability and Yawing 
Characteristics of an Untapered 45° Sweptback-Wing Model of 
Aspect Ratio 2.61. NACA TN 1668, 1948. 

9. Moller, Egon: Sechskomponenten-Messungen and Rechteckflli- 

geln mit V-Form und Pfeilform in einem grossen Schiebewin- 
kelbereich. Luftfahrtforschung, Bd. 18, Lfg. 7, July 19, 1941, 
pp. 243-252. 

10. Shortal, Joseph A.: Effect of Tip Shape and Dihedral on Lateral- 

Stability Characteristics. ‘ NACA Rep. 548, 1935. 

11. Bamber, M. J., and House, R. O.: Wind-Tunnel Investigation of 

Effect of Yaw on Lateral-Stability Characteristics. I — Four 
N.A.C.A. 23012 Wings of Various Plan Forms with and without 
Dihedral. NACA TN 703, 1939. 

12. Hollingworth, Thomas A.: Investigation of Effect of Sideslip on 

Lateral Stability Characteristics. Ill — Rectangular Low Wing 
on Circular Fuselage with Variations in Vertical-Tail Area and 
Fuselage Length with and without Horizontal Tail Surface. 
NACA ARR L5C13a, 1945. 

13. Irving, H. B., Batson, A. S., and Warsap, J. H.: Model Experi- 

ments on the Rolling Moment due to Sideslip of Tapered Wing 
Monoplanes. R. & M. No. 2019, British A.R.C., 1939. 

14. Multhopp, H.: Die Berechnung der Auftriebsvertcilung von Trag- 

flugeln. Luftfahrtforschung, Bd. 15, Lfg. 4, April 6, 1938, pp. 
153-169. (Available as R.T.P. Translation No. 2392, British 
M.A.P.) 


This formula is given in reference 14 along with other perti- 
nent formulas from which the lift coefficient, induced drag, 
and induced yawing moment may be determined for known 
loading ordinates. Reference 16 gives a formula for the 
lateral center of pressure in terms of the loading ordinates 
for m—M— 7. 


15. Crout, Prescott D.: A Short Method for Evaluating Determi- 

nants and Solving Systems of Linear Equations with Real or 
Complex Coefficients. Trans. AIEE, vol. 60, 1941, pp. 1235- 
1240. 

16. Multhopp, H.: The Application of the Wing Theory to Problems 

of Flight Mechanics. Reps, and Translations No. 860, British 
M.A.P. Volkendrode, Dec. 15, 1946. 



Table I 

CALCULATION SHEET FOR ASYMMETRICAL WING LOADS BY SEVEN-POINT WEISSINGER METHOD 

[A. X ; A ; tan A ] 


© 

© 

© 

© 

© 

© 

© 

© 

'© 

09 

® 

@ 

© 


. © 

© 

© 

© 


V 

3? 

y—y 

risa- 

1 

1 s=“ 

1 

v+v 


(y+y)- 

D 

ar v 

y tan A 

© + ® 

© 

<M 

© 


2 tan A 

v © 

X @ 

tan A 
X® 

®-l-@ 

@2 


V® 4 +® 

v®+© 

0. 9239 

1. 0000 

-0. 0761 

- 13. 1372 

0. 5198 

0. 0058 

3. 7013 

0. 8536 











. 9239 

. 9239 

0 

CO 

. 5412 

0 

3. 4142 

. 8536 







0 



CO 

. 9239 

. 7071 

. 2168 

4. 6132 

. 6131 

. 0470 

2. 6601 

. 8536 











. 9239 

. 3827 

. 5412 

1. 8477 

. 7654 

. 2929 

1. 7071 

. 8536 











. 9239 

0 

. 9239 

1. 0824 

1. 0824 

. 8536 

. 8536 

. 8536 











. 7071 

1. 0000 

2929 

-3. 4143 

. 5858 

. 0858 

2. 9142 

. 5000 










■ 

. 7071 

, 9239 

2168 

-4. 6132 

. 6131. 

. 0470 

2. 6001 

. 5000 











. 7071 

. 7071 

0 

CO 

. 7071 

0 

2. 0000 

. 50p0 







0 



00 

. 7071 

. 3827 

. 3244 

3. 0823 

. 9176 

. 1053 

1. 1876 

, 5000 











. 7071 

0 

. 7071 

1. 4142 

1. 4142 

. 5000 

. 5000 

. 5000 











. 3827 

1. 0000 

6173 

-1. 6199 

. 7232 

. 3811 

1. 9118 

. 1464 











. 3827 

. 9239 

5412 

-1. 8477 

. 7654 

. 2929 

1, 7071 

. 1464 











. 3827 

. 7071 

3244 

-3. 0823 

.9176 

. 1053 

1. 1870 

. 1464 











. 3827 

. 3827 

0 


I. 3086 

0 

. 5858 

. 1464 







0 



CO 

. 3827 

0 

. 3827 



. 1464 

. 1464 

. 1464 












■ 1 3[1— (1— X)y] 

ar~ A(l+X) 

>’ Computation unnecessary as It will load to indotorinlnant value. 
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Table I — Concluded 

CALCULATION SHEET FOR ASYMMETRICAL WING lT>ADS BY SEVEN-POINT WEISSINGER METHOD 


CO 

CO 

o 





@ 

© 

® 

□ 




B 


© 

® 

@ 



^31 

®x® 

X© 


® , 
© 

®x® 

X® 

V 

n 


/n.M 

©+■©+© 

La(v,m) 

V 

n 

Qr.n 

(0.0625) 

® 

®x® 

c,,„ 


v®+@ 






1 

i 

0 



( d ) 

■ ( d ) 

( d ) 

( d ) 

( d ) 

( d) 

( d ) 

( b ) 





1 

i 

1 

1.0000 


1 

i 




10. 4524 

■ 






T 

i 

2 

2. 0000 


1 

2 




-3. 6954 







1 

i 

3 

-I. 0000 


1 

3 




0 







1 

i 

4 

.4142 


( d ) 

( d ) 

( d ) 

( d ) 

( d ) 

( d ) 

(-) 






2 

2 

0 

1. 0000 


( d ) 

( d ) 

( d ) 

( d ) 

( d ) 

( d ) 

(■*) 






2 

2 

1 

. —2. 8284 


2 

1 




-2. 0000 


( b ) 





2 

2 

2 

0 


2 

2 




5. 6568 







■a 

2 

3 

2. 8284 


2 

3 




-2. 0000 







■ 

D 

4 

-1.0000 


( d ) 

( d ) 

( d ) 

( d ) 

( d ) 

( d ) 

c) 







i ' 

0 

-. 4142 


( d ) ■ 

( d ) 

( d ) 

( d ) 

( d ) 

' ( d ) 

c d ) 






3 

3 

1 1 

1. 0000 


3 





0 







ga 

3 

gSj 

-2. 0000 


3 

j| 




- 1. 5308 


( b ) 





| 

3 


-1. 0000 


3 

SB 

I ■ 



, 4. 3296 








3 

IB 

2. 4142 


( d ) 

("> 

( d ) 

( d ) 

(■') 

( d ) 

,( d ) 


b Computation unnecessary as it will lead to indeterminant value. 
■ For column @=tan A. 

4 No value exists. 
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Table II 


COLUMNS FOR INSERTION IN TABLE I WHEN MAKING 15-POINT WEISSINGER CALCULATIONS 

FOR ASYMMETRICAL WING LOADS * 






0. 98079 


. 98079 


. 98079 


. 98079 


. 98079 


. 98079 


. 98079 


. 98079 


. 98079 


. 92388 


. 92388 


. 92388 


. 92388 


. 92388 


. 92388 


. 92388 


. 92388 


. 92388 


. 83 


. 83147 


. 83147 


. 83147 


. 83147 


. 8 


. 83147 


. 83147 


. 83147 


. 70711 


. 70711 


.70711 


. 70711 


. 70711 


. 70711 


.70711 


1. 00000 


. 98079 


. 92388 


. 83147 


.70711 


. 55557 


. 38268 


. 19509 


0 


1. 00000 


. 98079 


. 92388 


. 83147 


. 70711 


. 55557 


. 38268 


. 19509 


. 98079 


• 92388 


. 83147 


. 70711 


. 38268 


. 19509 


0 


1. 00000 


. 98079 


. 92388 


. 83147 


. 70711 


. 55557 


. 38268 




2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4 

4 




-5. 02752 


2. 41408 


3. 53088 


n 

(0.03125) 

© 


( b ) • 


( b ) 


41. 0060 



*■ Columns ® to © of table I should be recalculated using columns © and © given here. Calculations corresponding to ij=5 should 
be deleted in columns <8, and @. 
e No value exists. 
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Table II — Concluded 

COLUMNS FOR INSERTION IN TABLE I WHEN MAZING 15-POINT WEISSINGER CALCULATIONS' 

FOR ASYMMETRICAL WING LOADS * 


® 

® 

© 

© 


© 

© 

© 


® 

V 

? 

V 

71 

A* 


V 

n 

(0.03125) 

® 

c r t n 

0. 70711 

0. 19509 

4 

4 

7 

2. 16480 

4 

7 


-0. 3170 

. 70711 

0 

4 

4 

8 

-1. 00000 

( b ) 

( b ) 

( b ) 

00 

. 55557 

1. 00000 

5 

5 

0 

-. 66816 

00 

( b ) 

00 

00 

. 55557 

. 98079 

5 

5 

1 

1. 41408 

5 

1 


0 

. 55557 

. 92388 

5 

5 

2 

-1. 69568 

5 

2 


3308 

. 55557 

. 83147 

5 

5 

3 

2. 41408 

5 

3 


0 

. 55557 

. 70711- 

5 

5 

4 

-4. 82848 

5 

4 


-3. 7936 

. 55557 

. 55557 

5 

5 

5 

-. 41408 

5 

5 


9. 6214 

. 55557 

. 38268 

5 

5 

6 

5. 69568 

5 

6 


-3. 7324 

. 55557 

. 19509 

5 

5 

7 

-3. 41440 

5 

7 


0 

. 55557 

0 

5 

5 

8 

1. 49664 

C b ) 

( b ) 

00 

00 

. 38268 

1. 00000 

6 

6 

0 

. 41440 

( b ) 

( b ) 

00 

■ 00 

. 38268 

. 98079 

6 

6 

1 

86720 

6 

‘ 1 


0550 

. 38268 

. 92388 

6 

6 

2 

i. 00000 

6 

2 

• 

0 

. 38268 

. 83147 

6 

6 

3 

-1. 29760 

6 

3 


-. 2976 

. 38268 

. 70711 

6 

6 

4 ' 

2. 00000 

6 

4 


0 

.38268 

. 55557 

6 

6 

5 

-4.35936 

6 

5 


-3. 3590 

. 38268 

. 38268 

6 

6 

6 

-1. 00000 

6 

6 


8. 6590 

. 38268 

. 19509 

6 

6 

7 

6. 52416 

6 

7 


-3. 1164 

. 38268 

0 

6 

6 

8 

-2.41408 

00 

( b ) 

( b ) ' 

OO 

. 19509 

1. 00000 

7 

7 

0 

-. 19904 

00 

( b ) 

00 

CO 

. 19509 

. 98079 

7 

7 

1 

. 41408 

7 

1 


0 

. 19509 

. 92388 

7 

7 

2 

-. 46912 

7 

2 


0518 

. 19509 

. 83147 

7 

7 

3 

. 58592 

7 

3 


0 

. 19509 

. 70711 

7 

7 

4 

-. 82848 

7 

4 


2286 

. 19509 

. 55557 

7 

7 

5 

1. 41440 

7 

5 


0 

. 19509 

. 38268 

7 

7 

6 

-3. 53088 

7 

6 


-2. 9356 

. 19509 

. 19509 

7 

7 

7 

-2. 41408 

7 

7 


8. 1572 

. 19509 

0 

7 

7 

8 

5. 02752 

( b ) 

( b ) 

00 

00 


» See footnote a,p. 13. 
b No value exists. 
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Table III 


SAMPLE CALCULATION OF ROLLING SPAN LOADING BY 7-POINT WEISSINGER METHOD ■ 


1 


1 



H 

@ 

® 

m 

. H 

1 


(°) 

V 

n 


fn.it 

@+®+© 

LaO,a) 

V 

71 


(0.0625) 

® 

®x© 



1 

i 

0 

-2.4142 

2. 1413 

w 

( d ) 

( d ) 

( d ) 

( d ) 

( d ) 

( d ) 

1 

i 

i 

L 0000 

2. 4107 

i 

l 

0. 9429 

0. 3049 

0. 2875 

10. 4524 

10. 7399 

1 

i 

2 

2. 0000 

2. 6406 

i 

2 

. 2559 

. 3049 

. 0780 

-3.6954 

-3.6174 

1 

i 

3 

-1. 0000 

2. 7037 

i 

3 

. 0913 

.3049 

. 0278 

0 

.0278 

1 

i 

4 

.4142 

2. 7141 

( d ) 

( d ) 

C d > 

C d ) 

C d ) 

C d ) 

( d ) 

2 

2 

0 

1. 0000 

1. 2620 

( d ) 

« 

( d ) 

C d ) 

( d ) 

C d ) 

( d ) 

2 

2 

1 

-2.8284 

1. 5532 

2 

1 

L 7852 

. 2538 

.4531 

-2.0000 

-1. 5469 

2 

2 

2 

0 

2. 4060 

2 

2 

1. 6620 

. 2538 

. 4218 

5. 6568 

6. 0786 

2 

2 

MM 

2. 8284 

2. 6356 

2 

3 

. 0064 

. 2538 

.0016 

-2.0000 

-1. 9984 

2 

2 

4 

-1.0000 

2. 6606 

( d ) 

( d ) 

(.*) 

( d ) 

t d ) 

( d ) 

C d ) 

3 

3 

0 

-. 4142 

.9232 

( d ) . 

(“) 

( d ) 

( d ) 


C d ) 

( d ) 

3 

3 


1. 0000 

1. 0075 

3 

1 

. 2899 

.2029 

.0588 

0 

.0588 

3 

3 

2 

-2. 0000 

L 4234 

3 

2 

2. 2828 

.2029 

.4632 

-1. 5308 

-1.0676 

3 

3 

3 

-1. 0000 

2. 3780 

3 

3 

1. 4761 

.2029 

.2995 

4. 3296 

4. 6291 

3 

IB 

4 

2. 4142 

2. 5167 

( d ) 

( d ) 

( d ) 

C d ) 

( d ) 

C d ) 

( d ) 


* Calculations o£ columns © to @ of table I are believed, to need no explanation and so are omitted. Calculations are carried out here from column @ through the 
determination of the loading coefficients for a wing of aspect ratio 3.5, sweep of 30°, and taper ratio of 0.3. 
a Sample calculations for column© 

When columns® and @=1, column @=(-2.4142) (2.1413)+(1.<X100) (2.4107)+ (2.0000) (2.&406)+(-1.0QQ0) (2.7037)+(0.4142) (2.7141) =0.0429 
When column @=1 and column ®=2, column @=(1.0000) (2.1413) +(-2.8284) (2.4107)+(0) (2.5M6)+(2.S2S41 (2.7037)+C-1.0000) (2.7141) =0.2559 
When column @=3 and column ®=2, column® =(1.0000) (0.9232) +(—2.8284) (1.0075)+(Q) (1.4234) +(2.8254) (2.3780)+(— 1.0000) (2.5167) =2.2828 
« Equations obtained from column @ 

10.7399Gi-3.6174Gi+0.0278G>=0.9239 
-1.5469G 1 +6.07S6Gi-1.9984G 3 =0.7071 
0.0588Gr-1.0676Gj+4 6291Gs=0.3817 
Gj= 0.1518 
Gj-0.1964 
Gi— 0.1260 
a No value exists. 

































































































































































